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ABSTRACT 

The spectra of A r affine Toda field theories with imaginary coupling constant, are investi- 
gated. Soliton solutions are found, which, despite the non-unitary form of the Lagrangian, 
have real classical masses and are stable to small perturbations. The quantum corrections 
to the soliton masses are determined, to lowest order in %. The solitons have the same 
spectrum as the fundamental Toda particles; a feature that is preserved in the quantum 
theory. 
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Introduction 



Amongst the conformal field theories in two dimensions, there exists a much simpler 
class of theories, the so-called 'minimal theories', whose spectra consist of a finite number 
of representations of the - possibly extended - conformal algebra [1] . For example, there 
are the theories where the conformal algebra is extended to the W r+ \ -algebras [2]. For 
each r there exists a whole set of such theories whose central charges are conveniently 
parameterized in terms of a constant a: 



For c > r - so a is real - the theories are unitary, but not minimal. For c < r - so a 
is purely imaginary - the situation is more complicated: when a 2 is rational the theory 
is minimal, since a finite set of representations decouple from the complete Hilbert space. 
Furthermore, when a 2 takes the particular rational values 



where k is a positive integer, the theory is, in addition, unitary. So, although a generic 
theory for which c < r is non-unitary, there exists a set of measure zero unitary theories, 
when the coupling constant takes the rational values of (2). 

This paper constitutes a study of the analogue of this phenomenon in the space of 
massive two-dimensional theories. In fact, we shall consider the massive theories which are 
thought to describe the W r +i -conformal field theories perturbed by the particular relevant 
operator which preserves the integrability of the theory [3]. With a as in (2), the operator 
has dimension {k + l)/(r + k + 1). A deformation which preserves integrability is very 
special, and suggests that the resulting massive theory may be tractable. 

There is evidence to suggest that these integrable theories are actually affine Toda field 
theories [4]. The case when r = 1, which is the sine-Gordon theory, has been studied in 
some detail [5], however, the more theories for more general r have not been the subject of 
much study - one of the reasons being that they are manifestly non-unitary. We show that 
although the theories are non-unitarity, there does exist a sector of soliton excitations which 
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appears to be unitary. This suggests that the soliton spectrum is the massive remanent 
of the primary conformal fields of the minimal conformal field theory which describes the 
ultra-violet scaling limit of the theory. 

The Lagrangian density of an affine Toda field theory - for a simply-laced Lie algebra 
- can be written in a number of ways: for instance 

C = \ ((cW - (d^) 2 ) - ^ J> (e^ - 1) . (3) 

7 j=o 

The field <fi(x,t) is an r = rank(g) vector, where g is a simply-laced Lie algebra, and 
consequently an algebra of type A, D or E. The inner products are taken with respect to 
the Killing form of g restricted to the Cartan subalgebra. The ck/s, for j = 1, . . . , r are the 
simple roots of g; a is the extended root (minus the highest root of g). The fact that the 
extended root is included in the sum, distinguishes the affine theories from the non-affine 
ones. The ct/s are linearly dependent: 

r 

^^•0^ = 0, n = 1, (4) 

3=0 

which serves to define the integers ny, 7 and m are constants. The quantum field theory 
defined by (3), is then thought to describe the integrable perturbation of the VF((7)-algebra 
conformal field theory 1 , with the identification 7 = \f^na. Notice that the particular value 
of 7 required to reproduce the minimal theories is actually imaginary. This implies that 
the Lagrangian density is no longer real - except for g = Ai, when it is the sine- Gordon 
theory. Prima facie, this would imply that the quantum field theory defined by (3) would 
be non-unitary. To mirror the situation at the critical point, we propose that when 7 2 /47r 
takes the rational values of (2), the Toda theory admits a unitary restriction. With this 
as the ultimate goal we now pursue the more modest aim of understanding the spectrum 
of the Toda theory with imaginary coupling constant 7. 

Before proceeding it is worth briefly reviewing the properties of the quantum theory for 
real values of 7. The approach that has been adopted in the literature entails considering 
(3) in the weak coupling limit, and then comparing conjectured forms for the full S- 
matrix with the weak coupling perturbation expansion [6]. The reason why this strategy 
is feasible is because the theory is expected to be integrable also at the quantum level, and 

1 The WV+i-algebra theories are those for which g = A r . 
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so the S'-matrix should be factorizable and also, in this instance, purely elastic [6]. The 
weak coupling expansion involves expanding (3) in powers of 7 about the minimum of the 
potential at <j> = 0; the first two terms are 

V(cf>) = ^ - ^) 2 + ^ E n 3 {a d .<pf + 0( 7 2 ). (5) 

j=o j=0 

In this limit, the theory describes a set of r particles - the fundamental Toda particles - 
whose masses are given by 2m\/C, where the £'s are the the eigenvalues of the matrix 

r 

J2 n 3 a 3® a 3> ( 6 ) 
3=0 

interacting via a set of three-point couplings defined by the term of order 7. The higher 
order terms describe corrections at various powers of 7. The first remarkable fact is that 
in the one-loop approximation all the masses are renormalized by an overall constant 
[6]. (This result is peculiar to the simply-laced theories.) On the basis of this, the first 
conjecture is that the mass ratios are preserved to all orders in 7. Furthermore, there are 
no additional states, so the quantum theory preserves the classical spectrum. For example, 
in the A r theories the classical masses are 

m a = 2m sin (j^p^j » a =1,2, ...,r, (7) 

and in the quantum theory the bare mass m is replaced by a renormalized mass niQ. 
Recently, the conjectured S'-matrix of the fundamental particles has been written in an 
elegant universal way using properties of the root system of the Lie algebra [7] . 

If 7 is imaginary, say 7 = i(3 for real /?, then a naive analytic continuation of the 
particle S'-matrix would lead one to conclude that the spectrum discussed above cannot 
be complete because of the appearance of new poles on the physical strip, indicating the 
existence of new particles states. Moreover, some the the residues of the S'-matrix, at the 
poles on the physical-strip, change sign - the hallmark of non-unitarity. In fact, we shall 
argue that the spectrum of the theory is actually much richer than this, because with 
imaginary coupling constant the theory also admits solitons. The solitons are present at 
the classical level, and correspond to static solutions of the equations of motion, in some 
frame of reference, with finite energy. It is straightforward to see why solitons will exist: 
with imaginary coupling constant the potential of the theory is periodic 4> ~ <fi + (2tt//?)w, 
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where uj is an element of the weight lattice A* of g. So the field <fi takes values on the 
maximal torus of g. Solitons correspond to field configurations which wind around the 
non-trivial cycles of the maximal torus, as x goes from — oo to oo. Thus the solitons are 
'kink' solutions, generalizing the situation for the sine-Gordon theory. What is not clear 
- since the potential is complex - is that the classical solitons have a real positive energy, 
and that they are stable to small perturbations. Remarkably both the former and latter 
are true for the A r theories, as we shall show. It turns out that the mass spectrum of 
solitons is equal to (7), up to an overall factor, although the solitons carry an additional 
quantum number which measures the topological charge. We further show that quantum 
corrections do not destroy this structure. 

The paper is organized as follows. In §1, the classical theory is considered, multiple 
soliton solutions are constructed and the classical soliton masses are calculated. In §2 the 
quantum corrections to the soliton masses are calculated. §3 contains some comments. 

1. The Classical Solitons 

In this section we will construct classical soliton solutions of the A r theories described 
by (3) and calculate their masses. 

1.1 Multiple Soliton Solutions 

Our approach to constructing the soliton solutions follows the methods developed 
by R. Hirota for finding explicit expressions for the multi-soliton solutions of integrable 
equations [8]. There are alternative formalisms for finding these solutions, most notably 
the inverse scattering method, however, for this particular task the Hirota method is far 
superior; its elegance being manifested in a subtle kind of non-linear superposition principle 
for soliton solutions, which easily allows one to construct multiple soliton solutions. The 
method is intimately connected with the tau-function approach to integrable theories, of 
which an up to date survey may be found in [9]. 

Before more general theories are treated, it is useful to recall some properties of the 
sine-Gordon theory, which is obtained by setting r = 1. With the definition ip = y/2(j), the 
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equation of motion, which follows from (3), is 

(9 t 2 -^)V+^-sin^ = 0. 

It has been known for a long time (see [10] for example) that this equation admits a soliton 
solution of the form 

i> = | tan" 1 ( e "(*-«*-0) , (1.1) 

where a, i> and £ are constants, and cr 2 (l — t> 2 ) = 4m 2 . The anti-soliton corresponds to 
the negative of the above. (The direct method was applied to the theory by R. Hirota in 

[11]-) 

We now turn to the more general theories. The crucial step of the Hirota method 
is a change of variables which converts the equation of motion to an equation of 'Hirota 
bilinear type' (to use the parlance of the subject) . The appropriate change of variables is 
not obvious; fortunately we are guided by the change of variables which is used for the 
Toda lattice equations, which are identical to the Toda field theory equations except that 
in the former the field only depends on t [12]. We are led to consider 

1 r 

P j=0 

where, for the moment, we are considering a general simply-laced theory. The new variables 
are the tau- functions for this theory. In terms of the r/s, the equations of motion become 

f f. T . _f 2 - t" t + t /2 \ 

I>( ' 1+ ' U'Ey-II^'-' (1-3) 

j=0 \ 3 J j=0 k=0 

(Here, dot denotes a t-derivative and prime an x-derivative.) This equation can be decou- 
pled to give the collection of equations 

r 



T jTj -TJ- TjTj + T? = m\ J] - r? ) • (1.4) 



The matrix with elements Ajk is the adjacency matrix 2 for the Dynkin diagram of the 
affine algebra In decoupling (1.3), the coefficient of r 2 has been fixed without loss of 
generality. 



2 This is the matrix which encodes the connectivity of the Dynkin diagram, such that Ajj = 0, 
and Ajk = 1 if the j th and k th spot are joined, i.e. aj - at = — 1. 
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When g = A r the equations are 

fjTj - if - t"t + rf = m 2 (T j _ 1 T j+1 - r?). (1.5) 

The label on Tj is to be understood modulo r + 1, owing to the periodicity of the extended 
A r Dynkin diagram. Those familiar with the Hirota method will identify (1.5) as being 
of a 'bilinear type', whereas, on the contrary, for the other simply-laced algebras this does 
not seem to be the case. We do not consider these other algebras further, however, we will 
make some comments about the -D4 theory in §3. 

To find multi-soliton solutions of the equations, the Hirota method starts by consid- 
ering a series expansion in some arbitrary parameter e, which will be set to 1 at the end 
of the calculation. So 

r J =l + er J (1) +e 2 rf ) + ---. (1.6) 

A solution describing N solitons results from terminating (1.6) by setting rj a ' > = for 
a > N. The fact that the series can be terminated is a crucial property of Hirota bilinear 
equations. Furthermore, one of the nice features of this approach is that the solutions have 
a very universal form for all integrable theories. The one-soliton solution has rj"^ = for 
a > 1, and rj satisfies the linear differential equation 

fW- r W» = ^( r W +r <i) 1 -2r«), (1.7) 



and the condition 



which are solved by 



#WiW )a > (1.8) 



T V = exp{<T(x-vt-Z)}6 j , (1.9) 
for constants a, v, £ and Sj. The constants Sj have to satisfy the recursion relations 

S^S; 1 = 6,5-^, (1.10) 

and the periodicity requirement Sj+ r +i = Sj. It is straightforward to solve for Sj, yielding 

Sj=u j , (1.11) 

where u is an (r + l) th root of unity. So there are r non-trivial solutions for which u = 
u a = exp 2nia/ (r + 1). In addition the constants a and v, and the integer a G {1,2,..., r}, 
have to satisfy 



a 2 (l-v 2 ) =4m 2 sin / ( -— ^ ) = m 2 a7 (1.12) 



where the m a 's are the classical masses of the fundamental particles defined in (7). (1.12) 
is just an expression of the fact that the theory is Lorentz invariant; the faster a soliton 
goes the narrower it becomes. 

The one-soliton solution is therefore 

0(a) (M) = ~ E «i lQ S ( X + ^-^^l) ■ (1-13) 
P i=0 

The solution is in general complex, however, the asymptotic values of 4>^ as a; — > ±oo are 
real. If we assume that a > 0, as a; — > — oo — > 0, and as x — > oo 

27T 

0^— (A a modA), (1.14) 

where A a is the a th fundamental weight 3 and A is the root lattice of g. We define the 
topological charge a configuration to be 

t = — I dxd x (j)= — lim — lim 4>(x). (1-15) 

So the one-soliton solution cf>^ has a topological charge which is equal to the a th funda- 
mental weight, up to a root. The actual value depends on a choice for the branch-cuts of 
the logarithms in (1.13). A careful analysis of the possibilities shows that the possible topo- 
logical charges of the one-soliton solution (j)^ fill out the weights of the a th fundamental 
representation. 

When g = A\ the solution can be chosen to real by choosing exp <r£ = —i exp with 
£' real, i.e. 

Tl =r * = 1 - ie°( x - vt -t'\ (1.16) 

reproducing the sine-Gordon soliton solution (1.1). 

Multi-soliton solutions result from terminating (1.6) at a higher power of e. It might be 
thought that the complexity of the solution increases rapidly with the number of solitons; 
however, this is not the case due to a subtle form of non-linear superposition principle 
which underlies the Hirota method. The two soliton solution illustrates the situation. The 
solution depends on a set of data {a p , v p , £ p , a p }, with p = 1,2; v p are the velocities, £ p are 
the position offsets and a p are the size parameters of the solitons. It is useful to define 

= a p (x - v p t - £ P ) + log (uQ , (1.17) 



The fundamental weights are dual to the roots: otj ■ Xk = 5jk- 
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then the form of the two soliton solution, after setting e = 1, is 



^(1) /r.(2) , (1) , i (2) 

Tj = l + e*5 + + Ae*i . (1.18) 

Substituting this ansatz into (1.5) fixes the 'interaction constant' 

(<7l - (T 2 ) 2 - ((TiUi - CT 2 V 2 ) 2 ~ 4m 2 sin 2 ( -^r(ai - a 2 ) ) 

A = -j (. (1.19) 

(ai + cr 2 ) 2 - ( cri^i + a 2 v 2 ) 2 - Am 2 sin ( ^-(ai + a 2 )J 

The expression for the multi-soliton solution builds on this elegant form: for example, the 
three soliton solution is 



t j = 1 + e i + e + e j + A v ; e j j + 



A^e^+^f + A^e*™**™ + A™ A*- 2 ® e*™ +*? )+ *?\ 



(1.20) 



The general N soliton solution can be written in the compact form 



(N N \ 

IX + E /W^ , (i-2i) 
P=l P,9=l / 



W}=C 

where = logA( p<? ). The sum is over the 2 N possibilities for which ji p = or 1, for 
each p. 

The multi-soliton solution have a natural physical interpretation. They represent the 
histories of a set of solitons which scatter off each other. To make this more precise consider 
the two soliton solution. Suppose that £1 < £ 2 , v\ > v 2 , and o\ and cr 2 are both positive. 
Focus on the solution in the vicinity of the first soliton, i.e. x ~ v\t + £1. Initially, that is 
in the limit t — > — 00, the solution is approximately 

t ~ 1 + e ^ , 

and finally, z.e. as t — > 00, the solution is approximately 

rj ~ e i ( 1 + Ae i J . 

In both the limits, the solution represents an isolated soliton, the only difference is that 
the final 'position offset' has been displaced: £1 1 — > £1 + In A. The fact that the velocities 
and form are unchanged by the collision is a characteristic feature of soliton scattering in 
an integrable field theory, reflecting the existence of the infinite set of integrals of motion. 
So asymptotically the solution represents two isolated solitons. What is not so apparent is 
that the topological charges of the solitons can be altered by the scattering (although the 
representations, i.e. a\ and a 2 , are fixed). To see this requires a rather detailed analysis, 
that will not be pursued here, of the positions of the branch-cuts of the logarithms. 
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1.2 Classical Soliton Masses 



The purpose of this section is to calculate the classical masses of the solitons that were 
constructed in the last section. The result is highly significant because it transpires that 
the ratios of the soliton masses are equal to the ratios of the masses of the fundamental 
particles (7). The result is that the one-soliton solution of (1.13), with a topological charge 
t = X a mod A, has classical mass 



M?=4^sin(-^)=^„. (i22) 



To find the masses of the soliton solutions, either the energy 

poo I 2 r \ 

M 1 {v) = J dx l^ + ^~Eni(e^-l) j , (1-23) 
or the momentum 

/oo 
dx4>-4>\ (1.24) 
-oo 

of a soliton must be calculated. (As usual, in the above 7(f) = (1 — v 2 )~ 1 ^ 2 .) The 
momentum is somewhat easier to calculate. In doing so it is useful to introduce the 
weights of the (r + 1) -dimensional representation 4 of A r , and to define <f>j = e,- • <fi. The 
momentum is then 

/oo r+1 n+ 1 „<^. + 

-°° j=i j=i j<t>j- 

where <f>j± are the asymptotic values of <j>j(x, t) as x — > ±00. Using the explicit expression 
for the one-soliton solution 0( o )(x, t) in (1.13) one finds 



(fij = — — ^- In 



i(3 dt \T~j-i 



(1.26) 



The weights of the (r + l)-dimensional representation are ej, for j = 1 to r + 1, with 

-\r- 
■>3 = 



ej = 0; in terms of which the roots are ctj = ej — e J+ i, with eo = e r +i- 
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where 5 = u a and $ = a(x — vt — £). This expression can be rewritten in terms of 4>j using 
the expressions 

l + *e» 1-.5-1 ^ 1+Jie* / 1 — <5 -1 10 

and 

Using these (1.26) becomes 

The utility of this expression is that <j>j has been written solely in terms of 4>j, and hence the 
integrals in (1.25) can now be performed. The constant terms in (1.29) do not contribute 
because X^=i 4>j = 0- Therefore 



j=i V / 



"3+ 



(1.30) 



The evaluation of this expression is straightforward using the asymptotic expressions for 



1 x — > —00. 



The following masses are found 



7(u)u tC^ 2 /3 2 1 



where the relation (1.12) was employed to eliminate a. 



2. Quantum Corrections to the Soliton Masses 



In this section we calculate the first quantum corrections to the soliton masses. We 
employ the semi-classical WKB method described in [13] for the kinks, or solitons, of the 
(p 4 and sine- Gordon theories - see also the review in [10]. 
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The dimensionless expansion parameter in our theory is h/3 2 /m 2 , and so the expan- 
sion in h coincides with the weak coupling expansion. For the sine-Gordon theory it is 
remarkable that the semi-classical WKB method, which gives the quantum corrections to 
the first order in h/3 2 /m 2 , is exact! The reason why the WKB method is exact for the 
sine- Gordon theory is, presumably, a consequence of its integrability. It seems plausible 
that the same is true for the affine Toda theories. In any case, we compute the mass cor- 
rections following exactly the same method as was used in [13] for the <f> and sine-Gordon 
theories - the latter being included in our calculation as the special case for which r = 1. 

From now on we set h = 1. The WKB method of [13] is straightforward to apply in 
our situation since the one-soliton solution is, in its rest frame, time-independent. The idea 
is to compute the zero-point energy of the small oscillations around the classical solution. 
The sum over all the modes can be done when an infra-red regulator is introduced; this is 
most easily achieved by putting the theory in a box. The zero-point energy of the vacuum 
must then be subtracted. The resulting expression is independent of the size of the box, L, 
as L — > oo; however, the final expression has an ultra-violet divergence. This divergence is 
not a problem of the soliton solution per se, rather, it is just a manifestation of the fact that 
the bare mass needs to be renormalized. This is achieved simply by normal-ordering the 
Lagrangian. The extra correction removes the divergence and the resulting finite residue 
which remains is then the mass correction. 

To begin with, we must analyse the small perturbations around the one-soliton solu- 
tion. If we write 4> = + rj, where 0( a ) is the static soliton solution in (1.13)- with v = 
and a = m a - and expand in powers of rj, to first order in rj we find 

r 

(d 2 - d 2 ) 77 + m 2 e^^ajiaj ■ rj) = 0. (2.1) 

j=0 

For a generic field theory, it is impossible to find the solutions to the linearized equation 
of motion around a classical solution of the equations of motion. However, we are in 
the fortunate position of studying an integrable theory for which we can find the modes 
explicitly. Asymptotically, as x — > ±oo, (2.1) becomes 

r 

(^ - d 2 x ) rj + m 2 otj (a 3 - - V ) = 0. (2.2) 

3=0 

This is nothing but the linearized classical equation of motion of the theory around the 
'vacuum' field configuration = 0. Since rj is an r-dimensional vector there are r indepen- 
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dent solutions of the full equation (2.1), with asymptotic form 



Vb (x,t) — > e^+^+^^wfo, 6 = l,-..,r, (2.3) 

with frequency 

U 2 = k 2 +ml, (24) 

which correspond to the fundamental particles modes of the Toda theory. We now proceed 
to find the exact solutions to the linearized equation (2.1). With (2.3) in mind, we are led 
to consider the two-soliton solution in eq. (1.18). The first soliton will be the static one- 
soliton solution that we are expanding around, the second will act as a small perturbation 
of the first. We take = ik and o^v 2 = —iv and treat exp o" 2 ^2 as a small parameter, in 
which case 

1 - s 



1 x -v St. 

j=o 

where r, corresponds to the one-soliton solution and 



^ = e*r+^ 1)+ *r. (2.5) 



The explicit expression for the 6 th mode is therefore 

m(k;x,t) = -Tg j> ( V/Tl — ) ^ n ™ J , (2-6) 



1 + A^e™^)\ j Akx+vt) 



+ u J ae m a (x-0 



with frequency z/ 2 = k 2 + m 2 . The interaction parameter A was found in (1.19), and for 
the present solution it is 

A , k) = ml + m 2 b -m 2 a _ h -2im a k 
m 2 a + m 2 b -m 2 a+b + 2im a k' 

It is straightforward to verify that the solutions (2.6) then have the correct asymptotic 
behaviour as predicted in (2.3). 

In addition to the solutions found above, there is a zero-mode solution (having zero 
frequency) of the form 

(,s) 
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This solution exists because the theory is translationally invariant and so there exists a 
whole family of one-soliton solutions which only differ in the position of the centre-of- 
mass. For the WKB method, it turns out that one can ignore this mode [13]; a more 
sophisticated time-dependent treatment shows that its effect is simply to ensure that the 
quantum soliton behaves as a relativistic particle, i.e. has energy levels E = \fM 2 + P 2 , 
where M is the quantum mass of the soliton and P is its momentum. 

Before we proceed to calculate the quantum mass of the soliton, it is worth pointing 
out that since there are no modes with imaginary frequency the solitons are classical stable. 
This is quite surprising given the fact that the potential of the theory is complex. 

In order to calculate the zero-point energy of the modes it is necessary to introduce 
some finite boundary conditions, whence the modes become discrete and it is possible to 
enumerate them. The appropriate boundary conditions in this case are rj(x = 0) = rj(x = 
L) = 0, where the centre-of-mass of the soliton lies between inside the box, and L is 
much larger than the size of the soliton ~ m~ l . The solutions satisfying these boundary 
conditions are 

r) b (k m ;x,t) - rj b (-k m ;x,t), (2.9) 

with 

k m L + Pb(k m ) = 7rm, m G Z. (2-10) 

where Pb(k) = Im(log A(k)). Notice that Pb(—k) = —pb(k). 

The quantum correction to the soliton mass is then obtained by taking the zero-point 
energy of the modes around the soliton solution and subtracting the zero-point energy of 
modes around the vacuum, which satisfy (2.10) with p = 0. That is 

AM « = ^EE (\A- + < ~ \l^ m + Pb(k m )/Lf + ml\ , (2.11) 

6=1 m >0 ^ ' 

where k m satisfies (2.10). Since L will eventually be taken to infinity, we can expand in 
L _1 , the leading term in the sum being 

kpb(k) __ pb{k) ^ (2 . 12) 



a//c 2 + ml dk 

with tb(k) = a/A; 2 + m\. Now, we take L — > oo in which case we can replace the sum over 
m by an integral over k: 

£=i{f ^+0(i- 1 )V (2.13) 



m>0 
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Therefore, the mass correction is 

dk n sdeb(k) 



AM a = -W°° d ± Pb (k) 

6=1 



dk 



=4t{/»(*M*)ir. + £*e.w^} 



(2.14) 



Now as k — > ±00, eb(/c) — > and 

~" 2^fc ^ + 2m ' " m '" b " m ' +b ^ ' (2 ' 15) 
and thus the first contribution from (2.14) is 

r _. r 

5>6(A0e6(ACoo = — E ( 2m « + - " ™«+*) • ( 2 - 16 ) 
6=1 ma b=i 



The sum is straightforward to perform, the result being a contribution 

m a {r + 1) 
2tt ' 



(2.17) 



to the mass. As it stands, the second contribution in (2.14) is not well defined because 
the integral has a logarithmic divergence. The problem is not due to any special nature 
of the soliton solution, indeed it is to be expected and is due to the fact that we have not 
renormalized the Lagrangian (3). We now pause to consider this aspect in more detail. 

Mercifully, the renormalization process in a two-dimensional field theory is straight- 
forward. The divergences can simply be removed by normal-ordering the Lagrangian. 
Working to lowest order in j3 2 and introducing an ultra-violet cut-off A 



. jfiar* := e W<*i-4> (l + ^£ (C a • aj ) 2 A a + 0(/3 4 ) ) , (2.18) 

where 



0=1 



f A dk 1 , 

A ° = Lr«-^wT^: (2 ' 19) 

and ( a is the a th eigenvector of X^=o a j ® a i °^ eigenvalue (m a /m) 2 . Hence, to lowest 
order in /3 2 , the renormalized Lagrangian is equal to 

L = \ {{d^f - (d x <P) 2 ) + ± ]T(m 2 + dm 2 ) (e* a '-+ - l) , (2.20) 

f j=o 
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with 

r 

dm] = \{m(3f ]T(Ca • a/Ao, (2.21) 

a=l 

and we have added a constant to bring the energy of the vacuum to zero. The presence 
of the additional terms clearly changes the energy of the soliton solution. These changes 
must be added to the quantum correction to the soliton mass. It might be thought that we 
should now consider the one-soliton solution for the renormalized Lagrangian; fortunately, 
this is unnecessary since the correction due to change in the form of the solution does not 
occur at the lowest order in /3 2 , precisely because the soliton solution is an extremum of 
the action. 

Taking into account the counter-term from the renormalization, the final expression 
for the mass correction for the a th soliton is 



AM a = - m - (n + 1) + y [ A %{ k ) dpb{k) 

y.TT < J I » Air 



J=0 

We now pause to calculate the integral 



2tt J _ a 4-7T dk 

6=1 J A 

-2 / dx (e l/3aj '^ - l) dm]. 



(2.22) 



/oo 
-oo 



dx fe i/3a '-*c«) - l) . (2.23) 



Firstly 

r 

e i0<*i-* w = Y[ r ~ ayah = Tj ~ 1 p +1 , (2.24) 
where Tj = 1 + uj J a e rna( ^ x ~^ is the static one-soliton solution. Using (1.5) we see that 

I _ T " T -+T /2 
TTT, 2 t 2 

, « ( 2 - 25 ) 

2 W J e m a (x-£) 



(1 + uj J ae m ^ x -^) 2 
The integral (2.23) may now be calculated explicitly yielding 



/m a \ 2 1 



•^f. (2.26) 

m z 
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The important point about the resulting expression for the integral is that it is in- 
dependent of j. The contribution to the mass correction (2.22) from the renormalization 
counter-term can now be simplified: 



i r r°° / - \ 

~7vH dx (^ arha) - A dm 



= \ m * Yl J2^ b ■ a i) 2 ^ ( 2 - 27 ) 

6=1 j=0 



i f m b\ 2 f dk 1 



where we used the fact that is an eigenvector of Y^j=o a j ® a j °^ eigenvalue (m^/m) 2 . 
The expression for the mass correction is now 



AM,- ma{r + 1) 



1 a 



-A r A dk j i fmb\ 2 1_ 



(2.28) 



V + ^( (m 2 +m 2_ m 2 +b)2+4m 2 fc2) j' 

Using the fact that 

r 

^(2m^ + 2m 2 - ™ 2 a _ b - m 2 a+b ) = 2m 2 a (r + 1), (2.29) 



^0 =Tl E Q i® a i J =X;«i-«i = 2(r + l), (2.30) 



6=1 

and 

/rat y 
V m / 

6=1 \j=0 / j=0 

one can easily verify that for large \k\ the integrand behaves as k~ 2 and so, as required, 
the logarithmic divergence is cancelled by the renormalization counter-term. We can now 
let the ultra-violet cut-off A tend to infinity. 

The remaining contribution is given in terms of a convergent integral: 

AM a = m a + la] , (2.31) 



where 



2tt 



dk ^ J [To , „_ 2 4(m^ + m 2 . - m 2 +{) ) 



OO 



777- 



1 



(2.32) 
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Unfortunately, it is not possible to evaluate the integral I a analytically; except when r = 1 
in which case the integrand is zero. Actually, this case is of interest since it corresponds 
to the sine-Gordon theory, and we find 

M q = M cl + AM 

_8m_ V2m (2-33) 

So if we define M cl (/3 2 ) to be the classical mass of the soliton, then to lowest order in the 
quantum corrections 

M q = M cl (/3' 2 ), (2.34) 

with 

f = T^W.< < 2 ' 36 > 

the well-known result for the sine-Gordon quantum soliton mass. Different quantization 
approaches show that this result is exact [14]. 

Returning to the more general theories, we can solve the integral numerically. For a 
large number of values of r we have found that I a is actually independent of a. The first 
10 values are presented in the table 5 . 



r 


I 


1 





2 


-1.814 


3 


-1.571 


4 


-2.672 


5 


-2.721 


6 


-3.620 


7 


-3.792 


8 


-4.593 


9 


-5.577 



Table: the integral I(r) 



The important point is not the actual value for I a , but the fact that it is indepen- 
dent of a. This means that to lowest order the soliton spectrum just receives an overall 



5 Numerical integration was performed using Mathematica 
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renormalization and the masses are still proportional to the masses m a . Our final results 
is 

M2 = m a + /) + 0(/3 2 ), (2.36) 

where f3' is defined in (2.35). I is the constant which only depends on r. 

3. Discussion 

We have demonstrated that Toda theories with imaginary coupling constant admit 
soliton solutions, whose topological charges lie in the set of weights of the fundamental 
representations of the algebra. Surprisingly the solitons have real classical masses and 
are stable to perturbation. The gross spectrum of the solitons is the same as that of the 
fundamental Toda particles. This feature is not spoiled by the quantum theory, to lowest 
order in h, and with the sine-Gordon example in mind, we can speculate that, as a result 
of the integrability of the theory, the masses in (2.36) are also exact. 

In the sine-Gordon theory, the soliton and anti-soliton can form bound states - the so- 
called breathers. The masses of the breathers can be found using more sophisticated time- 
dependent WKB methods [13]. Remarkably the ground state of the breather spectrum is 
identified as the fundamental particle of the sine-Gordon theory. It is clearly of interest 
to repeat the calculation of ref. [13] and establish whether the more general Toda theories 
admit breather states. One would also like to determine the soliton S-matrix. Since the 
theory is integrable we expect this S-matrix to be factorizable, however, the solitons carry 
topological charge and consequently we do not expect the S-matrix to be elastic, i.e. we 
expect non-trivial scattering in the space of topological charge. It appears that the correct 
S'-matrix has the form of a product of the minimal S'-matrix, which determines the pole 
structure, and a quantum group .R-matrix which acts in the space of the topological charge; 
details of this will be presented elsewhere. When the S-matrix has been determined, it 
should then be possible to determine out whether the solitons can be decoupled from the 
non-unitary sector of the theory. 

We have only considered the A r affine Toda field theories and it would be interesting 
to know whether the construction of the solitons can be extended to the other simply- 
laced theories. Initial considerations for the affine theory based on D4 indicate that soliton 
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solutions do exist, and are also associated to the fundamental representations of the algebra. 
Furthermore, the classical masses of the solitons are consistent with the conjecture 

Oh 

Mf = -^m a , (3.1) 

where h is the Coxeter number of the algebra, and the m a 's are the masses of the fun- 
damental particles of the theory; that is 2my / C, where ( is one of the r eigenvalues of 
(6). 

Soliton-type solutions of the A 2 affine Toda field theory are also considered in [15]; 
however, they are not the same as the solutions considered here, since they do not satisfy 
the classical equations of motion. 

I would like to thank Mark de Groot for some helpful comments on the manuscript. 
This work was supported by NSF PHY90-21984. 
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